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Abstract 

We developed a new method for full kinematical reconstruction of the tt system near 
its threshold at future linear e~^e~ colliders. In the core of the method lies likelihood 
5^ I fitting which is designed to improve measurement accuracies of the kinematical variables 

that specify the final states resulting from tt decays. The improvement is demonstrated 
by applying this method to a Monte-Carlo tt sample generated with various experimental 
effects including beamstrahlung, finite acceptance and resolution of the detector system, 
etc. In most cases the fit brings a broad non-Gaussian distribution of a given kinematical 
variable to a nearly Gaussian shape, thereby justifying phenomenological analyses based 
on simple Gaussian smearing of parton-level momenta. The standard deviations of the 
resultant distributions of various kinematical variables are given in order to facilitate such 
phenomenological analyses. A possible application of the kinematical fitting method and 
its expected impact are also discussed. 
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1 Introduction 



The discovery of the top quark ilj at Tevatron has completed the standard-model (SM) list 
of matter fermions. In spite of its subsequent studies thereat, our knowledge on its proper- 
ties is still far below the level we reached for the lighter matter fermions. A next-generation 
e+e~ linear collider such as JLC j2], having a facet as a top-quark factory, is expected to 
allow us to measure top quark's properties with unprecedented precision, thereby improving 
this situation dramatically. Such precision measurements may shed light on the electroweak- 
symmetry-breaking mechanism or hint beyond-the-SM physics or both. 

Being aware of the opportunities provided by the linear collider, a number of authors 
have so far performed interesting analyses on the measurements of top quark properties [SI 
miSlEI- They can be classified into two categories, i.e., those near the tt threshold, mainly 
focused on physics contained in the threshold enhancement factor, and those in open-top region, 
searching for anomalies in production and decay vertices, both of which play important roles 
and complement each other. 

In those analyses, feasibility studies on form factor measurements have been done mainly 
in the open-top region [HIHlini- In the meantime, it has been conceived that, in view of the 
energy upgrading scenario of the e^e~ linear collider, measurements of top form factors in the 
tt threshold region are also important. The top quark physics is expected to commence in the 
threshold region at the early stage of the collider operation and full exploration of the machine 
potential in that phase is crucial for the project design. It has, therefore, been repeatedly 
stressed that a realistic simulation study is in desperate need to clarify feasibility of precision 
measurements of form factors at the tt threshold. Besides, form factor measurements in the tt 
threshold region have some favorable features: availability of well-controlled highly polarized 
top sample ^01 no need for transformation to t or t rest frames because both t and i 
are nearly at rest; as far as the decay form factor measurements of an on-shell top quark are 
concerned, the center-of-mass energy does not matter. 

In order to thoroughly carry out such analyses for real data, we need a sophisticated 
method to kinematically reconstruct events as efficiently and as precisely as possible. This 
is, however, highly non-trivial in practice, due to finite detector resolutions, possible missing 
neutrinos in the final states, and various background contributions. Furthermore, care has to be 
taken when imposing a kinematical constraint on the masses of the t and t quarks because they 
cannot be simultaneously on-shell below the threshold. We thus need to further explore the 
potential of the e+e~ linear collider and extend the past studies [3 |Hl 1121 US] to the threshold 
region, in order to make maximum use of the linear collider's advantages: clean experimental 
environment, well-defined initial state, availability of highly polarized electron beam, possibility 
of full parton-level reconstruction of final states, etc. 

In this paper, we thus aim at developing an efficient method for full kinematical recon- 
struction of the ti system near its threshold in e~^e~ annihilation and clarifying the accuracy 
to which various observable will be measured. We develop a likelihood fitting method which 
is especially designed to improve measurement accuracies of kinematical variables of the par- 
ticles originating from the tt sample in the threshold region. Moreover, some of the analysis 
techniques developed here are expected to be useful for the analyses in the open-top region. 

Our study should also provide important information to the current line of phenomeno- 
logical studies on top quark physics at linear colliders. In fact, there have been a number of 
theoretical studies on measurements of the top-quark production and decay form factors using 
the e+e~ ti process El El El- However, many of these analyses assumed either the most 
optimistic case or the most conservative case with respect to the kinematical reconstruction of 
event profiles. In the former case, one assumes that the momenta of all the particles (including 
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t and W) can be determined precisely, while in the latter case, one uses only partial kinemat- 
ical information, e.g. the direction of b, and the energy and momentum of i. In this work we 
will provide realistic values of resolutions with which individual kinematical variables can be 
measured. 

In Sec. 2 we briefly review our simulation framework. Sec. 3 is devoted to top quark 
reconstruction in the lepton-plus-4-jet mode, where two subsections recapitulate basic strategy 
and procedure, respectively. In Sec. 4 we explain our kinematical reconstruction using a likeli- 
hood fitting method. Then we discuss a possible application of this method and its expected 
impact in Sec. 5. Finally, Sec. 6 summarizes our results and concludes this paper. 

2 Framework of Analysis 

For Monte-Carlo-simulation studies of ti productions and decays, we developed an event gen- 
erator that is now included in physsim-2001a ^1], where the amplitude calculation and phase 
space integration are performed with HELAS/BASES [T^ ITB] and parton 4-momenta of an event 
are generated by SPRING [TO]. In the amplitude calculation, initial state radiation (ISR) as well 
as S- and P-wave QCD corrections to the ti system [TTl HH] are taken into account. Parton 
showering and hadronization are carried out using JETSET 7 . 4 ^H] with final-state tau leptons 
treated by TAUOLA [201 in order to handle their polarizations properly. 

In this study, the top-quark (pole) mass is assumed to be 175 GeV and the nominal center- 
of-mass energy is set at 2 GeV-above the 15* resonance of the ti bound states. This energy is 
known to be suitable for measurements of various properties of the ti system at threshold jT2]. 
We will assume an electron-beam polarization of 80% in what follows. Effects of natural beam- 
energy spread and beamstrahlung are taken into account according to the prescription given in 
|12j . where the details of the beam parameters are also described. We have assumed no crossing 
angle between the electron and the positron beams and ignored the transverse component of 
the initial state radiation. Consequently, the ti system in our Monte-Carlo sample has no 
transverse momentum. Under these conditions we expect 40k ti events for 100/6~^. 

The generated Monte-Carlo ti events were passed to a detector simulator (JSF Quick 
Simulator j2l]) which incorporates the ACFA-JLC study parameters (see Table. The 
quick simulator created vertex-detector hits, smeared charged-track parameters in the central 
tracker with parameter correlation properly taken into account, and simulated calorimeter 
signals as from individual segments, thereby allowing realistic simulation of cluster overlapping. 
It should also be noted that track-cluster matching was performed to achieve the best energy- 
flow measurements. 



Detector 


Performance 


Coverage 


Vertex detector 


= 7.0 © (20.0/p) / sin^/2 q 




cos 6 


< 0.90 


Central drift chamber 


ap^/pT = 1.1 X 10-V © 0.1 % 




cos 9 


< 0.95 


EM calorimeter 


(Xe/E = 15 % /VE © 1 % 




cos 9 


< 0.90 


Hadron calorimeter 


as/E = 40 % /VE © 2 % 




cos 9 


< 0.90 



Table 1: ACFA study parameters of the JLC detector, where p, pt, and E are measured in 
units of GeV. 
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3 Event Selection 



3.1 Basic Reconstruction Strategy 

Since the top quark decays almost 100% into a b quark and a W boson, the signature of 
a tt production is two b jets and two W bosons in the final state. These W bosons decay 
subsequently either leptonically into a lepton plus a neutrino or hadronically into two jets. 
According to how the W bosons decay, therefore, there will be three modes of final states: (1) 
six jets, where both of the W^s decay hadronically, (2) one lepton plus four jets, where one 
of the W^s decays leptonically and the other hadronically, and (3) two leptons plus two jets, 
where both of the W^s decay leptonically. 

In order to reconstruct the momentum vector of the top quark, we will use the lepton- 
plus-4-jet mode, for which we can reconstruct the t(f)-quark momentum as the momentum sum 
of the b{b) jet and the two jets from the hadronically-decayed W~^{W~), while we can tell the 
charge of the hadronically-decayed W from the charge of the lepton. 

In the lepton-plus-4-jet mode, two of the four jets are b(b) jets directly from the t{t) 
quarks, while the other two are from the W boson that decayed hadronically. Therefore, if one 
can identify the b and b jets, remaining two jets can be uniquely assigned as decay products 
of the W boson. The other W boson can be reconstructed from the lepton and the neutrino 
indirectly detected as a missing momentum. Remaining task is then to decide which b{b) jet 
to attach to which 14^-boson candidate, in order to form t{t) quarks. Since the t{t) quarks 
are virtually at rest near the threshold, a b{b) jet and the corresponding W boson fiy in the 
opposite directions. We can thus choose the correct combination by requiring the b{b) jet and 
the W boson be approximately back-to-back. 

In reality, however, 6(6)-quark tagging is not perfect and can be performed only with 
some finite efficiency and purity: there could be more than two 6(6)-jet candidates in a single 
event. In addition, b and b quarks can be emitted in the same direction. In such a case, a 
wrong combination could accidentally satisfy the back-to-back condition. These facts sometimes 
prevent us from uniquely assigning each jet to its corresponding parton, resulting in multiple 
solutions for a single event. Moreover, the leptonically-decayed W is poorly reconstructed in 
practice, since the neutrino momentum is strongly affected by ISR, beamstrahlung, as well as 
other possible neutrinos emitted from the b{b) jets. In order to overcome these difficulties, we 
will need some sophisticated method. We defer discussion of such a method to the next section 
and examine here the extent to which the aforementioned basic reconstruction strategy works. 

3.2 Event Selection Procedure 

The lepton-plus-4-jet-mode selection started with demanding an energetic isolated lepton: Ei > 
18 GeV and -Ei40cone < 18 GeV, where Ei is the lepton's energy and -Ei40cone is the energy sum 
of particles within a cone with a half angle of 14° around the lepton direction excluding the 
lepton itself.* When such a lepton was found, the rest of the final-state particles was forced 
clustering to four jets, using the Durham clustering algorithm |22]- Two-jet invariant mass was 
then calculated for each of the six possible combinations and checked if it was between 65 GeV 
and 95 GeV, in order to select a jet pair which was consistent with that coming from a M^-boson 
decay. For such a jet pair the remaining two jets, at the same time, had to be identified as b(b) 
jets, using fiavor tagging based on the impact parameter method. The hatched histogram in 
Fig. ^ is the 2-jet invariant mass distribution of all the possible pairs out of the four jets, while 

*The Ei cut was chosen to be the kinematical limit for the lepton from the W ^ tv decay. On the other 
hand, the cone-energy cut was optimized to achieve high purity, while keeping reasonable efficiency. 
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the solid histogram being that with the 6-tagging. It is seen that this procedure dramatically 
improved the purity of the W boson sample. It should also be stressed that these selection 
criteria are very effective to suppress background processes such as e~^e~ — >■ W~^W~ and provide 
us with an essentially background-free tt event sample. 




Figure 1: Invariant mass distributions of the 2-jet systems reconstructed as W^-boson candi- 
dates. Hatched and solid histograms correspond to before and after double 6-tagging, respec- 
tively. The locations of the W mass cuts are indicated with arrows. 

The remaining task is to decide which b{b) jet to associate with which W candidate. For 
a 6(6)-jet candidate, the right W boson partner was selected by requiring the back-to-back 
condition as described above. Fig. El is a scatter plot of the acoplanarity angles of the two 
possible b-W systems where horizontal and vertical axes are the angles of h-Wty and b-W2-jet 
system, respectively. b-W pairs having 9acop{b-w) < 60° was regarded as daughters of the t{i) 
quarks. 
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Figure 2: Scatter plot of the acoplanarity angles corresponding to two b-W systems, where 
horizontal and vertical axes are angles of b-Wny and b-W2-jet systems, respectively. 

The selection efficiency after all of these cuts was found to be 15% including the branching 
fraction to the lepton-plus-4-jet mode of 29%. 
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4 Kinematical Fit 



The event selection described above yields a very clean ti sample. As noted above, however, 
the sample is still subject to combinatorial backgrounds, if we are to fully reconstruct the final 
state by assigning each jet to a corresponding decay daughter of the t or t quark. We thus need 
a well-defined criterion to select the best from possible multiple solutions. It is also desirable 
to improve the measurement accuracies of those kinematical variables which are suffering from 
effects of missing neutrinos (such variables include momenta of b, b or the neutrino from a W 
itself). 

The ti system produced via e~^e~ annihilation is a heavily constrained system: there are 
many mass constraints in addition to the usual 4-momentum conservation. At e~^e~ linear 
colliders, thanks to their well-defined initial state and the clean environment, we can make full 
use of these constraints and perform a kinematical fit to select the best solution and to improve 
the measurement accuracies of the kinematical variables of the final-state partons. 

4.1 Parameters, Constraints, and Likelihood Function 

For the lepton-plus-4-jet final state, there are 10 unknown parameters to be determined by the 
fit: the energies of four jets, the 4-momentum of the neutrino from the leptonically-decayed W 
boson, and the energies of the initial-state electron and positron, provided that the jet directions 
as output from the jet finder are accurate enough, the error in the 4-momentum measurement 
of the lepton from the leptonically-decayed W can be ignored, and that the transverse momenta 
of the initial-state electron and positron after beamstrahlung or initial-state radiation or both 
are either negligible or known from a low angle e~^/e~ detector system '''(see Fig. EI). 




Figure 3: Schematic diagram showing parameters and constraints relevant to the kinematical 
fit described in the text. The boxed parameters are unknown and to be determined by the fit. 

The requirements of 4-momentum conservation and the massless constraint for the neu- 
trino from the leptonically-decayed W reduce the number of free parameters to 5. We choose, 
as these free parameters, the energies of the four jets and the initial longitudinal momentum 
(the difference of the energies of the initial-state electron and positron). 

These five unknown parameters can be determined by maximizing the following likelihood 
function: 

L = iI[Pi {E~\ Ef)) . Pr^, ■ Pr^_ ■ Pr,, ■ P^s. (1) 
/^i 

^In addition, there will be some finite transverse momenta due to a finite crossing angle of the two beams. 
These transverse momenta are, however, known and can be easily incorporated into the fit. 
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where P^^ is a resolution function for jet / and is Gaussian for / = 1 and 2 (jets from the 
hadronically-decayed W) as given by the detector energy resolution. For / = 3 and 4 (jets 
from the b and b quarks) the resolution function is the same Gaussian convoluted with the 
missing energy spectrum due to possible neutrino emissions. For the two W bosons in the final 
state, we use a Breit-Wigner function Pp,^ instead of 5-function-like mass constraints. is a 
weight function coming from ISR and beamstrahlung effects. This distribution was calculated 
as a differential cross section as a function of the energies of initial-state electron and positron, 
taking into account the ti threshold correction as described in Sec. El 

The remaining factor, Pr^^-, controls the mass distribution of the t and i quarks and 
has been introduced to take into account the kinematical constraint that the t and t cannot 
be simultaneously on-shell below threshold (see Fig. |3] which shows Pr^j distribution below 
ti threshold). Pp^^ distribution is a dynamics- independent factor which is extracted from the 
theoretical formula for the threshold cross section. 
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Figure 4: Pp^- distribution below ti threshold. 



4.2 Results 

We performed the maximum likelihood fit for the selected sample. The maximum likelihood 
fit provided us with a well-defined clear-cut criterion to select the best solution, when there 
were multiple possible solutions for a single event: we should select the one with the highest 
likelihood. 

Figs. IHJ-a) and -b) are the reconstructed W mass distributions for the leptonically and 
hadronically-decayed W bosons, respectively, before (hatched) and after (solid) the kinematical 
fit. The figures demonstrate that the Breit-Wigner factors (Pp^^) in the likelihood function 
properly constrain the W masses as intended. 

Fig. El-a) plots the reconstructed mass for the t{i) decayed into 3 jets against that of the 
i{t) decayed into a lepton plus a 6 jet, before the kinematical fit. The strong negative correlation 
is due to the fact that the neutrino from the leptonically-decayed W is reconstructed as the total 
missing momentum. Figs. El-b) and -c) are the projections of Fig. El-a) to the horizontal and 
vertical axes, respectively, showing systematic shifts of the peak positions. Figs.El-d) through - 
f) are similar plots to Figs.|Hl-a) through -c) after the kinematical fitting, while Figs.EJ-g) through 

"i-This is in contrast with the result in |12| . where a quite tight set of cuts was imposed upon the reconstructed 
W and t masses, and consequently their peak shifts were less apparent at the cost of significant loss of usable 
events. The goal of this study is to establish an analysis procedure to restore those events which would have 
been lost, by relaxing the tight cuts while keeping reasonable accuracy for event reconstruction. 
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Figure 5: Reconstructed W mass distributions for (a) leptonically and (b) hadronically-decayed 
W bosons, before (hatched) and after (sohd) the kinematical fitting. Note that the vertical 
scale for the hatched area is enhanced by a factor of five for illustration purpose. 

-i) are corresponding distributions of generated values (Monte-Carlo truth). The kinematical 
fit sent most of the events to the L-shaped region indicated in Fig. l^-d), as it should, and made 
the distribution look like the generated distribution shown in Fig.|Hl-g). Consequently, the peak 
shifts observed in the Figs. IHl-b) and -c) have been corrected as seen in Figs. [HJ-e) and -f). There 
are, however, still some small fraction of events left along the minus 45° line. These events 
were so poorly measured that it was impossible to restore. The cut (angled region) indicated in 
Fig. ini-d) allowed us to remove them without introducing any strong bias on the reconstruction 
of the kinematical variables. 

Now the question is how the above constraints improve the parameters of the fit such as 
the energies of b and b jets, the direction and the magnitude of the missing neutrino from the 
leptonically-decayed W, on which we expect significant infiuences. Figs. [3-a) and -b) plot the 
difference between the reconstructed and the generated energies of the b (b) quark attached 
to the leptonically-decayed W and that of the b (6) attached to the hadronically-decayed W, 
respectively, before (hatched) and after (solid) the kinematical fit. The plots demonstrate that 
the kinematical constraints recover the energies carried away by neutrinos from the 6 or 6 decays. 
The kinematical fit brings broad non-Gaussian distributions into nearly Gaussian shapes. The 
standard deviations of the 6 or 6 jet energy distributions are approximately 3.5 GeV after the 
kinematical fit. 

The improvement is more dramatic for the direct neutrino from the leptonically-decayed 
W, which is reconstructed as the total missing momentum; see Figs. U^c) and -d) which show 
distributions of the difference of the reconstructed and generated neutrino energies (AEu) and 
directions {A6u). Again the kinematical fit makes the broad and skewed distribution of neutrino 
energies into a nearly Gaussian shape with a standard deviation of approximately 2.5 GeV. The 
fit also improves the angular resolution as shown in Fig. [3-d). The resultant angular resolution 
is ag = 2.9°, which was obtained by fitting Nq^ exp(— 6*^/2(7^) to the distribution. 

The improvements in these kinematical variables are refiected to the improvements in the 
reconstructed W energies and directions as shown in Fig. |H]-a) for the energy of the leptonically- 
decayed W, -b) for the hadronically-decayed W, and -c) for the direction of the leptonically 
decayed W. We can see dramatic improvements in all of these distributions, although the 
improvement in the direction of the hadronically-decayed W is less dramatic. Both the en- 
ergy resolution of the leptonically-decayed W and hadronically-decayed W are approximately 
2.4 GeV, the angular resolution of the leptonically-decayed W and hadronically-decayed W are 
2.4° and 1.7°, respectively. 
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Figure 6: Scatter plot of the reconstructed t{i) mass from 3 jets versus that from a lepton plus 
a 6-jet (a) before the kinematical fit, together with (b) its horizontal/fe^i/ and (c) vertical/3-jet 
projections, (d) through (f) are similar plots after the kinematical fit and (g) though (i) are 
corresponding plots for generated values. 
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Figure 7: Distributions of the difference of the reconstructed and generated energies of the b or 
6 jet attached to (a) leptonically-decayed and (b) hadronically-decayed W bosons, together with 
distribution of the difference of the reconstructed and generated (c) energies and (d) directions 
of the direct neutrino from the leptonically-decayed W, before (hatched) and after (sohd) the 
kinematical fit. 




Figure 8: Distributions of the difference of the reconstructed and generated energies of (a) 
leptonically-decayed and (b) hadronically-decayed W bosons, and (c) distribution of the dif- 
ference of the reconstructed and generated directions of the leptonically-decayed W, before 
(hatched) and after (solid) the kinematical fit. 
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Finally, we will examine the effects of the kinematical fit on the measurements of the 
direction and the magnitude of the top quark momentum. In Figs. IHl-a) and -b), the difference 
of the reconstructed and generated directions of the t or t quark is plotted against the generated 
top momentum, before and after the kinematical fit, respectively. We can see appreciable 
improvement by the fit. Nevertheless, since the top quark direction becomes more and more 
difficult to measure as the top quark momentum decreases, the resolution is still somewhat poor 
in the low momentum region. The angular resolution is largely determined by the reconstruction 
of the t or t decayed into 3 jets. Remember that the resolution improvements were less significant 
for the hadronically-decayed W, since the power of the constraints was used up mostly to recover 
the momentum information of the direct neutrino from the leptonically-decayed W and the 
energy resolution for jets from the W was left essentially unimproved. The improvement in the 
measurements of the top quark direction is mostly coming from the improvement in the 6 or 6 
jet measurement. By the same token, the effect of the fit on the measurement of the magnitude 
of the top quark momentum is also less dramatic compared to that on the leptonically-decayed 
W. The momentum and angular resolutions of the t or t quarks after the fit are approximately 
3.0 GeV and 5.5°, respectively. ^ 




Figure 9: The difference of the reconstructed and generated directions of the t or t quark 
plotted against the generated top momentum, (a) before and (b) after the kinematical fit. 

In the case of the 6-jet mode, for which there is no direct energetic neutrino from W^s, we 
can use the power of the constraints to improve the jet energy measurements. Consequently, 
we may expect more significant improvement in the top quark momentum measurement. 



5 A Possible Application 



We discuss a possible application of our kinematical reconstruction method. Let us consider 
measurements of the decay form factors of the top quark in the tt threshold region. We assume 
that deviations of the top-decay form factors from the tree-level SM values are small and 
consider the deviations only up to the first order, i.e. we neglect the terms quadratic in the 
anomalous form factors. Then the cross sections depend only on two form factors fi and 
in the limit rrib although the most general tbW coupling includes six independent form 
factors 



Wtb 



9W -I X 



7^ i'^PL 



-J2 -^R 



w 



t), 



(2) 



^ Since the distributions deviate from Gaussian shapes substantially off their peaks, these values should be 
taken as order of magnitude estimates. 
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where = (1 - 75)72 and Pr = (1 + 75)72. At tree level of the SM, f[ = 1 and = 0. 
A variation of fi changes only the normalization of the differential decay width of the top 
quark, whereas a variation of changes both the normalization and the shape of the decay 
distributions. Thus, we expect that the kinematical reconstruction is useful for disentanglement 
of the two form factors and in particular for the measurement of f^. For simplicity we assume 
f^ = l hereafter.^ Since transverse W (denoted as Wt) is more sensitive to than longitudinal 
W (Wl), our strategy is to extract Wt using the angular distribution oi W (in the rest frame 
of t) and the angular distribution of £ (in the rest frame of W). It is well known that Wt 
is enhanced in the backward region cos 6w ~ —1, where the angle 6w is measured from the 
direction of the top quark spin in the t rest frame. Also, we may enhance Wt by collecting 
£ emitted in the backward direction cos6'£ ~ —1, where the angle 6i is measured from the 
direction of —pt in the W rest frame. These features are demonstrated in Figs. ^1 We plot^' 
(a) the differential decay width for the decay of the top quark with a definite spin orientation 
dT{ti — >• Mu) /{d cos 9wd cos 9 e) for f^ = and (b) the difference of the differential widths for 
= 0.1 and for = 0. The plots show that we may measure f^, for instance, from the ratio 
of the numbers of events in the regions cos 9w, cos9i < and cos 9w, cos9i > 0. 




Figure 10: (a) Normalized differential decay width Af^"^ dV{t^ — >• Mp) / {d cos 9y/d cos 9 () for 
= 0. (b) Difference of the normalized differential decay widths for = 0.1 and for = 0. 
In both figures the differential widths are normalized by A/" = Ft x Br(iy £v) for = 0. 

In the last section, we showed the significant improvement of the reconstruction of the 
leptonically decayed W due to the kinematical fit. In order to see how the improvement affects 
the measurement of the distribution in question, namely that in Fig. ITUl-a). we compare the 
reconstructed and generated distributions before and after the kinematical fit, using the same 
Monte Carlo sample we used in the previous sections. Fig. ITTl-a) and -b) plot, for the selected ti 
sample, the reconstructed differential decay width normalized by the corresponding generator 
level distribution (a) before and (b) after the kinematical fit. It is clear from Fig. ITTl-a) that 
the measurement is biased towards high cos 6'^, which is because the energy of the leptonically 
decayed W tends to be overestimated so that the lepton from the W is often over-boosted. 
Fig. ITTl-b) demonstrates that the kinematical fit effectively removed such a measurement bias. 
We expect therefore that the kinematical fit will reduce possible systematic errors in the dif- 

^In order to determine fi simultaneously, we may, for instance, use independent information from the 
measurement of the top width fT^ . 

'I We used the helicity amplitudes given in I^Hl for calculating these differential decay widths. 
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ferential width measurement significantly, thereby improving sensitivity to f^.** 




Figure 11: Reconstructed differential decay width distributions corresponding to Fig. ITOl-a) 
normalized by the generator-level distribution (a) before and (b) after the kinematical fit. 

It is advantageous to investigate decay properties of the top quark in the tt threshold 
region as compared to the open-top region E ^ 2mt because of several reasons. First of all, 
the top quark can be polarized close to 100% in the threshold region 1^. which is a useful 
tool to sort out various form factors. This is clear in the above example. Furthermore, we 
are almost in the rest frame of the top quark. In the above example, the top quark is highly 
polarized in its rest frame. Hence, the event rate expressed in terms of cos6w and cos 6'^ is 
a direct measure of the amplitude-squared, \J2i=L,T-^(t'\ ~^ bWi) x A(Wi lu)\'^ (without 
phase-space Jacobian), which allows for simple physical interpretations of event shapes. We 
also note that we do not gain resolving power for the decay form factors by raising the cm. 
energy. This is in contrast with the measurements of the ti production form factors. 

6 Summary and Conclusions 

To make maximum use of future e^e~ linear colliders' experimental potential, the top quark 
reconstruction in the lepton-plus-4-jet mode has been studied under realistic experimental con- 
ditions of e^e~ —>■ tt process near its threshold. As a new technique to fully reconstruct ti 
final states, we have developed a kinematical fitting algorithm which aims to reconstruct the 
kinematical variables of top quarks and their offsprings more accurately. 

The missing energy carried away by neutrinos from bottom quark decays has been recov- 
ered by the kinematical fitting. However, the effects of the kinematical fitting on the top quark 
momentum are not as dramatic as we wanted. This is because the top quarks are almost at rest 
in the threshold region and therefore their momenta are difficult to measure. Moreover, in the 
lepton-plus-4-jet mode many constraints are used up by recovering the information on the neu- 
trino from leptonically-decayed W bosons. On the other hand, the remarkable improvements 
of the energy resolution of fe-jets and the angular and energy resolutions of leptonically-decayed 
W^s have been achieved by the kinematical fitting. These improvements should benefit the 
form factor measurements in general. As a possible application, we considered measurements 

**We can extract also from the distribution of £ energies measured in the laboratory frame without relying 
on the reconstruction of its parent W momentum. The sensitivity of the lepton energy distribution to /J^ is, 
however, estimated to be lower than that of the differential decay width. 
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of decay form factors including on which correct reconstruction of the leptonically-decayed 
W may have a large impact. 

As stated in Sec. ^ many theoretical studies on measurements of the top form factors 
assumed either the most optimistic case or the most conservative case with respect to the 
kinematical reconstruction of event profiles. Our analysis indicates that both assumptions 
are not realistic under actual experimental conditions. In this respect we emphasize that 
the kinematical fit brought often heavily skewed and broad distributions into nearly Gaussian 
shapes. Realistic phenomenological analyses using information of the decay particles from top 
quarks will then become possible by simply Gaussian-smearing parton-level momenta with the 
resolutions for the measurements obtained in this study. To be specific, the resolution for 
jet energy measurements is o"^. ~ 3.5 GeV after the kinematical fit for both the light quark 
jets from W boson decays and the bottom quark jets from t or t quarks. As for the energy 
resolution for the neutrino coming from the leptonically decayed W we have — 2.5 GeV. 
The energy resolutions for both of the leptonically and hadronically decayed W^s then become 

— 2.4 GeV, and the angular resolutions for the leptonically decayed W and the neutrino 
directly coming from it improve to 2.4° and 2.9°, respectively. Finally the momentum and 
angular resolutions for the the t or t quarks are approximately 3.0 GeV and 5.5°, respectively. 
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